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Abstract 

To describe the spin-orbital energy gap formation in the one-dimensional 
symmetrically coupled spin-orbital model, we propose a simple mean field 
theory based on an SU(4) constraint fermion representation of spins and or- 
bitals. A spin-orbital gapped phase is formed due to a marginally relevant 
spin-orbital valence bond pairing interaction. The energy gap of the spin and 
orbital excitations grows extremely slowly from the SU(4) symmetric point 
up to a maximum value and then decreases rapidly. By calculating the spin, 
orbital, and spin-orbital tensor static susceptibilities at zero temperature, we 
find a crossover from coherent to incoherent magnetic excitations as the spin- 
orbital coupling decreasing from large to small values. 
PACS Numbers: 75.10.Jm, 75.10.Lp, 75.40.Gb 
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It is a very interesting problem to look for exotic quantum magnetic states realized with 
the spin and orbital degrees of freedom. Since the discovery of new quasi-one-dimensional 
spin gapped materials Na 2 Ti 2 Sb 2 [1] and NaV 2 5 [2], there has been considerable interest 
in magnetic systems with orbital degeneracy [3-10]. It is believed that the unusual magnetic 
properties observed in these compounds can be explained by a simple two-band Hubbard 
model at quarter filling, and in the large Coulomb repulsion limit the effective Hamiltonian 
is simplified to a model of two symmetrically coupled spin -1/2 Heisenberg chains: [11-13] 

H — J y~] (Si • Sj+i + Tj • T i+ i) 

i 

+\/]T(S J -S m )(T l -T m ), (1) 

i 

where Sj and Tj denote the spin and orbital degrees of freedom, respectively. Here both 
coupling parameters J and V are assumed to be antiferromagnetic, and the model Hamil- 
tonian is SU(2)cg>SU(2) symmetric with an additional Z 2 symmetry in exchange between Sj 
and Tj. 

For this model, the critical point V = describes two independent isotropic Heisenberg 
spin-1/2 chains with gapless excitations. In the weak coupling regime V/J <C 1, it has been 
shown that the model describes a non-Haldane spin liquid where magnetic excitations are 
gapful but incoherent [14]. In the strong coupling regime [15] V/J ^> 1, however, a special 
point V/J = 4 has been identified where the Hamiltonian becomes SU(4) invariant [13], and 
it has been demonstrated by Bethe ansatz and effective field theory methods that the low- 
energy excitations are coherent, given by three branches of gapless elementary excitations 
[16,17]. However, as emphasized in Ref. [6], there is no renormalization group flow from the 
first to the second critical points due to the Zamolodchikov theorem on the central charge. 
In the Heisenberg limit (V = 0), the total central charge is c = 2, while c = 3 at the 
SU(4) symmetric point. This means a gapped phase is expected in-between the weak and 
strong coupling limits with gapless magnetic excitations. In particular, it has been shown 
[18] that when V/J = 4/3, the model has an exact ground state in which spin and orbital 
operators may form dimerized singlets in a staggered pattern, and such a matrix product 
state is doubly degenerate and gapped. However, so far it is not clear whether such staggered 
dimmerized singlets can represent the ground state in the whole gapped phase. On the other 
hand, a crossover transition from incoherent to coherent magnetic excitations is speculated 
as the spin-orbital coupling changing from small to large values within the gapped phase 
[6]. Since the perturbation treatment from either end can not provide a unified description, 
a non-perturbative interpolation scheme would be highly desirable. 

The purpose of the present paper is to develop a simple mean field (MF) theory based on 
the strong-coupling SU(4) symmetric limit and to describe such a spin-orbital gap formation 
along with the coherent-incoherent crossover of magnetic excitations. When V / J < 4 and is 
close to the strong-coupling SU(4) point, the spin and orbital static susceptibilities display 
very sharp coherent magnon peaks at the nesting wave vector 2k p = it/2, corresponding to 
a commensurate spin-density wave of period four lattice spacings. Away from the strong 
coupling symmetric point, a spin-orbital valence bond (VB) pairing interaction is present, 
leading to an energy gap in the quasiparticle excitations. The energy gap initially grows 
extremely slowly from V/J = 4 up to a maximum value near V / J ~ 1.265, and then decreases 
rapidly to a very small value near V/ J ~ 0.462. From the calculated spin and orbital static 
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susceptibilities, the coherent magnetic peaks around q = ±7r/2 are gradually suppressed 
and slightly shifted, and their spectral weights are transferred to the incoherent background 
around q = ±7r. Moreover, the present MF theory also provides the correlation spectra of 
the spin-orbital tensor operators, exhibiting further clear evidence of the nontrivial crossover 
from the weak to strong coupling limits of the model. 

First, up to a constant, the model Hamiltonian can be rewritten as: 



# = J c £ (2S, • S m + • T i+1 + 1) 

~Js X! (^Si ■ Si+i - 2) (^ 2 Ti • T i+1 - , 



(2) 



where the first part corresponds to an SU(4) spin-orbital symmetric model in the SU(4) 
fundamental representation [4-7,10], while the second part corresponds to a staggered SU(4) 
spin-orbital VB model [19-21,10] in which alternating sublattice sites transform according 
to the SU(4) fundamental and ant i- fundamental representations, respectively. Here the 
coupling parameters are regrouped into J c = (J/2 + V/8) and J s = (J/2 — V/8). 

In the strong coupling SU(4) symmetric point, V/J = 4, J c = J and J s = 0. Our 
MF theory will take this limit as a starting point, while the weak coupling limit V = 
corresponds to the case of J c — J s — J/2. In order to maintain the higher symmetry of the 
strong coupling limit and to characterize both spin and orbital degrees of freedom at the 
same time, an SU(4) constrained fermion representation is introduced, and its generators 
are given by Fp(i) = Cj^Ci^, satisfying the SU(4) Lie algebra 



(3) 



The four states we consider are |+,+ >,| — ,+ >, |+, — > and |— , — >, where the first index 
specifies the spin projection, while the second one is the orbital projection. It's then obvious 
that the spin and orbital operators are expressed in terms of these four-component fermions 

as: 



— C\ t iC it 2 + Cj j3 Ci t 4, 

Ft — c} tl c it3 + C\ 2 Ci^ 
Fi = Cj^Ci^i + Cj j4 C ii2 , 

rpz _ ^ (-1 , (-1 (-1 S~l] (-1 \ 

from which the following commutation relations can be proved 



(4) 



(5) 



Sf , Sj 

T + T~ 

1 ' 3 

qa rpft 
&i > I j 



2s;&. 

2T?5i 



± 



J i,3i 

0, a, (3 ■ 



rpZ T^i 

i ' 3 

x,y,z. 



(6) 
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It is thus demonstrated that the spin and orbital operators are two independent degrees of 
freedom and both of them obey their respective SU(2) Lie algebra. By imposing a local 
constraint Z)^ C/^Cj^ = 1, we can further prove that constraints S 2 = T 2 = 3/4 are 
satisfied, corresponding to the spin- 1/2 system with two-fold orbital degeneracy. Under the 
new representation, the model Hamiltonian is expressed as a quadratic form in terms of two 
composite operators 

H = - J c ]T : A\A t : -J s ]T B\B U (7) 

i i 

with 

Bi = [(Ci+l,4Q,l + Ci+l,lCi,4j 

_ (Ci+l,3Ci,2 + Ci+l,2Ci,3)] , (8) 

where describes a nearest neighbor VB hopping parameter, while Bi represents a nearest 
neighbor VB pairing parameter. The normal ordering has been chosen in the first term. 

To develop a MF theory, the nearest neighbor VB order parameters are defined by 
A c (i) = (Ai) and A s (i) = —(Bi). The model Hamiltonian is then decomposed into 

H = -J C E [&c(i)ClC i+1 ,, + H.c] +Xj2cla,, 

i,fi i,fi 

+j s J2 [a.(o (clc\ +lA + c\ A c\ +l:1 ) 

i 

-A s (i) (cJ 2 C/ +li3 + Ct 3 Cj +h2 ) + H.c' 



-XN + N 



J c \A c (t)\ 2 + J s \A s (i)\ 2 , (9) 



where a local chemical potential is first introduced to impose the local constraint and then 
it is replaced by a global value A keeping the translational symmetry. When the spatial 
uniformity of VB parameters are also assumed, in terms of a generalized Nambu spinor, 

*1 = (Cfc,l> Cfc, 2 > Cfc,3> Cl A , C_fc,l, C-fc,2, C , _ fej3 , C_fc i4 ) 

the MF model Hamiltonian can be rewritten in a compact form 

n = \ £ ^lH m /(A;)^ fc + AiV + N (j c A c 2 + J S A S 2 ) , 

k 

where H m/ (ife) = [A - A c (k)] Oi - A a (fc)n 2 , A c (ife) = 2J c A c cosk, A s (k) = (2J s A s sink) 
with Oi = cr 2 <E> o"o <H> o"o an d ^2 = &x ® o y <E> c^- The corresponding Lagrangian is given by 

L W = \ £ [iw n - H m/ (A;)] * fc (iw n ) + ... 

the Matsubara Green's function matrix is thus derived as 
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^ = ^%^ffl. ( 10 ) 

Then the fermionic excitation spectra with fourfold degeneracy are yielded 

e k = ±\j (A - 2 J C A C cos kf + (2 J S A S sin kf. (11) 

In the excitation spectra (with plus sign), the local energy minima appear at the specific 
momentum - the so-called Fermi momentum kp, where an energy gap opens up when A s ^ 0. 

Now consider the static properties at zero temperature. By filling in all states with 
negative energies, the ground state energy per site is evaluated as 

e g = - r — J(X- 2J c A c coskf + (2J S A S sin/c) 2 

J-7T 71 

+A + ( J C A C 2 + J S A S 2 ) . (12) 

By minimizing the ground state energy with respect to parameters A c , A s , and A, the saddle 
point equations are derived as: 



dk — (A — 2J C A C cos k) cos k 



-■ 

/7< 
-■ 



71 


>/(A 


- 2J c A c cos kf + (2J S A S 


sin kf 


dk 




2 J s sin 2 k 




71 




- 2J c A c coskf + (2J S A S 


sin kf 


dk 




(A - 2 J C A C cos k) 




7T 




- 2J c A c cos/c) 2 + (2J S A S 


sin kf 



A c , 
1, 

1. (13) 



In particular, when V/J — A, the self-consistent equations are easily solved, and we obtain 
A = 4J/V, A c = 2y/2/ V, and A s = 0. There are four degenerate gapless fermionic energy 
bands, different from the three bosonic elementary excitations obtained from the Bethe 
ansatz method [16]. However, as will be shown later, there are only three gapless collective 
(bosonic) excitations, so the physical conclusions are correct. The reason why four, instead of 
three, gapless modes show up in the fermion representation is similar to the weak coupling 
limit of the effective bosonization approach [6]. In that approach the charge excitation 
becomes gapped in the strong coupling limit due to an umklapp term, whereas the other 
three branches remain degenerate and gapless with marginally irrelevant interactions [6]. 
Away from the SU(4) symmetric point V/J < 4, numerical calculations can be performed 
and solutions to these self-consistent equations are derived: A and A c steadily decrease as 
the coupling parameter V / J is reduced, while A s gradually increases at the same time. 

The ground state energy per site is plotted as a function of the coupling parameter V / J 
in Fig.l. In order to compare two limiting cases, the corresponding ground state energies 
for A c 7^ 0, A s = and A c = 0, A s ^ are also plotted in the same figure as well. It 
has been found that the gapped phase (A c ^ and A s ^ 0) is smoothly connected with 
the strong coupling SU(4) symmetric gapless phase (A c ^ 0,A S = 0), and represents the 
possible lowest ground energy state in the parameter range of 0.462 < V/ J < 4. Our strong 
coupling MF theory is probably limited to this regime. 
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In Fig. 2 the fermionic quasiparticle spectra with both positive and negative energies are 
plotted in the range [— 7r/2,7r/2] for different couplings V/J = 3.85, 2.00, 1.265, 0.8. In 
particular, at a special value of (V/ J) c ~ 1.265, both spectra in the range [—hp, kp] become 
completely flat. Furthermore, for V/J > (V/J) c , the minimum point of the spectra appears 
at a finite Fermi momentum ±&f, while for V / J < (V/ J) c the minimum point of the spectra 
is shifted to zero. The special value of (V/J) c is very close to the exactly soluble point of 
V/J = 4/3 of the model Hamiltonian [18]. 

The Fermi momentum in the excitation spectrum mentioned above is given by 



kp = cos 1 



A(2J C A C 



(2J C A C ) 2 - (2J a A a y 



(14) 



which is also calculated as a function of V/J and plotted in Fig. 3a. It has been found that 
k F is almost fixed at n/A over a large range of 2 < V/J < 4, and then quickly decreases 
to zero near the point of V/J ~ 1.18. The energy gap opens up at momentum kp and is 
evaluated as 



Agap 



\ { 2JsAs)2 + pjij^j^H' (15) 



which is presented in Fig. 3b. In the range of 3 < V/J < 4, the energy gap is extremely 
small. This agrees with the recent density matrix renormalization group calculations showing 
exponentially slow gap opening [7], in contrast to the earlier results [3]. Only when V / J < 3, 
the energy gap starts to grow slowly up to a maximum near the critical value (V/J) c ~ 
1.265, and then decreases to a very small value near V/J ~ 0.462. The position of the 
maximum energy gap roughly corresponds to the condition of the dispersionless quasiparticle 
excitations, which is consistent with the analysis at the exactly soluble point [18]. 

In the present strong coupling MF theory, the spin-spin and orbital-orbital density cor- 
relation functions can simply be evaluated as well. The spin and orbital density operators 
Eq. (4) are re-expressed in terms of the generalized Nambu spinor 

s? = itftngtfi, t? = 

ft>S — a z ® CO ® fx, = CO ® CXq <E> <7y, 

= a z <S> CTq <S> &z, = a z ® °x ® Co, 

Q V T — a Q ® a y ® (J , Vt z T = a z ® a z ® (To- 
Then the spin and orbital density- density correlation functions are given by 

1 r dk 

X a x(q,iUm) = -jQp E J ^Tr[n a x G(k,iu; n )n°G(k + q,iu; m + iu; n )}, (16) 

where Q x = Qg for the spin and Q x = Qj. for the orbital. By inserting the Matsubara 
Green function, it is straightforward to prove the following relation, 

Xs(Q, IWm) = Xt(Q, iuJ m) = Xfa, iuj m), (17) 
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independent of the indices a = x,y, z. This shows that away from the SU(4) symmetry point, 
the spin and orbital rotational symmetry of SU(2)®SU(2) with an additional Z 2 symmetry 
in exchange between Sj and Tj is satisfied in the present strong coupling MF state. The 
resulting expression of x(<?, iw m ) is given by 



, m + iu n ) + (A - A c (k)) (A - A c (k + q)) + A s (k)A s (k + q) 



IjJn 

iu n (iu r , 



-e 



+ iu n ) -e{ +q 



(18) 



However, in the SU(4) spin-fermion representation, nine spin-orbital tensor operators 
associated with nonlinear collective excitation of both spin and orbital degrees of freedom, 
can be defined by = 2S?T?, which can also be expressed in terms of the generalized 
Nambu spinor as 



n x L x 

or 



C 2 



0~x 

(Ty 

&X 



X ■ 



nf 



(T ® (Ty 
(TO <8> (Tx 
(To® °~z [ 

(TO ® (Jy I 



yi 

(Ty, 



T Z . (19) 

Then it can be proved that the corresponding nine correlation functions (T T L?' /3 (r)L?' / V)) 
are the same, and equal to 

1 ^ r dk 

iu n (iu m + iu n ) + (A - A c (k)) (A - A c (k + q)) - A s (k)A s (k + q) 



(iuj n y 



(iu m + iuj n ) -e 



k+q 



(20) 



Compared to the spin and orbital density correlation functions, there is only a sign difference 
in front of A s (k)A s (k + q). In the conventional response function theory, the correlation 
spectrum Xhifli^m) represents a nonlinear collective excitations of both spins and orbitals. 
When V/J = 4, it is found that A s = 0, and then we have Xl (os^m) — X'sio.i^m) = 
XriQ-i^m)-! independent of their component indices of spin and orbital operators, which 
implies that the SU(4) symmetry is recovered. 

After summation over the Matsubara frequency and analytical continuation, we find that 
in both dynamic susceptibilities x(g, cj) and Xl(q,w), an energy gap exists at uj — 2A gap in 
the large parameter range of 0.462 < V/J < 4. As far as the relevant experiments [1,2] are 
concerned, however, the most important physical quantities measured are the corresponding 
static susceptibilities. In the zero frequency limit, the static susceptibilities x(o) an d Xl{q) 
are deduced to respectively 




/ 



+ 



dk f [n F (e k ) - n F (e k+q )] 
87r \ e k+q - e k 
[1 - n F (e k+q ) - n F (e k )] 

£k+q + £k 



1 + 



(A - A c (fc)) (A - A c (k + qj) ± A s (k)A s (k + q) 

e k^k+q 

(A - A c (fc)) (A - A c {k + q)) ± A s (k)A s (k + q) 

^k^k+q 
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We present the static spin and orbital susceptibility at zero temperature in Fig. 4 and the 
static spin-orbital tensor susceptibility at zero temperature in Fig. 5 for different coupling 
parameters V/J = 2.72, 2.00, 1.265, 0.80, and 0.50. 

At V/J — 2.72, since the energy gap in the elementary excitations is negligible, both 
static susceptibilities x(q) an d Xl(q) have sharp peaks at q — ±vr/2, corresponding to a 
commensurate spin- and orbital-density wave of period four. The period four arises because 
the spin and orbital chain are in their SU(4) fundamental representation and by "quadru- 
plicity" one needs four sites to form a singlet. Moreover, as the momentum goes to zero, 
x(q) approaches to zero, while Xl(q) to a constant. 

As V I J decreases further, a finite gap opens up, and the sharp peaks in both x(q) 
and Xl(q) spectra at q — ±7r/2 are strongly suppressed and broadened. Meanwhile, the 
peak positions are slightly shifted to lower momenta, indicating the possible presence of 
incommensurate density waves in the gapped phase. The spectral weights of the coherent 
quasiparticle peaks are gradually transferred to the incoherent excitations at q — ±7r for 
x(q) and at q = for Xl(q)- Thus, the static spin and orbital susceptibilities around q = 
±7r are enhanced, exhibiting a crossover from coherent to incoherent magnetic excitations. 
Some of these features have been speculated by an effective low-energy bosonization theory 
[6,14]. From our MF theory, the crossover is estimated to occur around the critical coupling 
V/J = 1.265. It seems to us that the present symmetric state is a good approximation 
of the genuine ground state to describe this crossover. Moreover, the spectrum of Xl(q) 
in Fig. 5 displays a further clear evidence of the crossover of the magnetic excitations from 
the strong to weak coupling limits of the system. In particular, when V/ J becomes smaller 
and smaller, a new peak structure clearly emerges at zero momentum. All these features of 
Xl(q) are new results obtained from the present non-perturbative theory. 

Finally, several remarks and comments are in order, i) The fermionic MF theory devel- 
oped in this paper mainly focuses on the energy gap formation of the spin-orbital coupled 
model, in particular, on the crossover of coherent-incoherent magnetic excitations in the 
gapped phase, ii) All the obtained results on the gapped phase are mostly consistent with 
the latest density matrix renormalization group calculation [7] and the effective bosonization 
theory [6]. One of the important results is that the energy gap opens up extremely slowly 
away from the SU(4) symmetric point, indicating that the phase transition from gapless 
to gapped phases may be of Kosterlitz-Thouless type. The gapped regime with coherent 
magnetic excitations shares some similarities with the Haldane gapped phase of the quan- 
tum antiferromagnetic spin-ladder model with a four-spin interaction [14]. So it is not clear 
whether the Lieb-Schultz-Mattis theorem is valid in the present spin-orbital coupled system 
or not. iii) The gapless phase at the SU(4) point is only characterized by a Fermi liquid be- 
havior in the present MF theory approximately, different from the Luttinger liquid behavior 
derived from the Bethe ansatz and effective field theories [16,17,6]. However, as far as the 
spin, orbital, and spin-orbital tensor collective excitations are concerned, all these physical 
results of the spin-orbital system — the correlation spectra are in good agreement with the 
exact solution. 

In conclusion, we have applied an SU(4) constraint fermion representation to a one- 
dimensional symmetrically coupled spin-orbital model, a spin-orbital gapped phase is gener- 
ated away from the strong coupling SU(4) symmetric point by a relevant spin-orbital pairing 
interaction. The energy gap of the elementary spin and orbital excitations grows extremely 
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slowly up to a maximum value and then decreases. The spin, orbital, and spin-orbital tensor 
static susceptibilities are also calculated at zero temperature, displaying a crossover from 
coherent to incoherent magnetic excitations as the spin-orbital coupling decreasing away 
from the SU(4) symmetry point. It is interesting to note that the present mean field theory 
provides a rather good description of the underlying physics over a large crossover region 
between the two critical points. It's likely that the SU(4) constraint fermion representation 
is appropriate for this system. 

We would like to thank Prof. Z. B. Su and X. Q. Wang for useful discussions, while L. 
Yu would like to thank A. Nersesyan for very helpful comments. This work is supported 
by NSF-China (Grant No. 10074036 and 10125418) and the Special Fund for Major State 
Basic Research Projects of China (G2000067107). 
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Figure Captions 

Fig. 1. The ground state energy per site as a function of the coupling parameter V/ J. For 
comparison, the corresponding ground state energies for A c ^ 0, A s = and A c = 0, A s ^ 
are also plotted by the dashed and dash-dotted lines, respectively. 

Fig. 2. The fermionic energy spectra are plotted in the range [— tt/2, +7t/2] for the 
coupling values: a) V/J = 3.85, 2.00 and b) V/J = 1.265, 0.80. 

Fig. 3. The Fermi momentum /cp and the corresponding quasiparticle gap as functions 
of coupling parameter V / J. 

Fig. 4. The static spin and orbital susceptibilities at zero temperature for different 
coupling values: a) V/J = 2.72, 2.00 and b) V/J = 1.265, 0.80, and 0.50. 

Fig. 5. The static susceptibility of the spin-orbital tensor at zero temperature for different 
coupling values: a) V/J = 2.72, 2.00 and b) V/J = 1.265, 0.80, and 0.50. 
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